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THE LIMITING VELOCITY IN FALLING FROM A GREAT HEIGHT. 

By Edwin Bid well Wilson. 

1. The fundamental characteristic of the vertical motion under gravity in a resisting 
medium is the approach to a final or limiting velocity U, whether the initial downward velocity 
is less or greater than U. The equations of motion are easily set up and integrated when ±he 
density of the medium is constant "and the resistance varies as the square of the velocity. 

A fact not so much stressed is the slowness of the approach to the limiting velocity U. 
For example the simple relation v 1 =2g7t, which neglects altogether the resistance of the atmos- 
phere, shows that the .height A of release from rest must he about 10,000 feet before the object 
will attain a velocity of as much as 800 ft./sec. even without the opposition of the air resistance; 
and when allowance is made for that resistance the height must he greater. If, therefore, 
terminal velocities of 900 or more feet per second, such as are customary with airplane bombs, 
are under consideration, it is only in" the case of fall from a great height (upward of 10,000 
feet) that the terminal velocity can be anywhere nearly approached. 

Here, however, another difficulty enters. The resistance of the air varies with the density 
of the medium and this variation should not be assumed a priori to be negligible in the case of 
fall from heights of upward of 2 miles. In the standard table of densities at different levels, 
the ratio p/p 0 of the density at altitude A to that at the earth has become about 0.74 at 1Q,000 
feet, 0.62 at 15,000, and 0.44 at 25,000. The question thus arises as to whether the changing 
resistance of the air may be taken into account in some satisfactorily simple way for the dis- 
cussion of vertical fall. 

It is the purpose of this report to give that simple treatment of the problem which I have 
in the past offered to my classes at the Massachusetts Institute of Technology; there are 
undoubtedly other solutions, perhaps equally simple or simpler, but I have seen none. 

2. The equation of motion is 

W Tt W+wty* 

if W be the mass; v the velocity upward, and epv 2 be the resisting force in pounds. With the 
substitution u=v > the equation reduces to it= — 2g and is linear, but with a variable 

coefficient p. The general solution may be indicated in the form 

ue ~ W** = -2g j 'e ~ TrJ" ^ dh+C 

In order to perform the integration it is necessary (unless graphical methods are used) to 
have an expression for p as a function of 7i that is sufficiently accurate and at the same time trf 
such form as to make integrable the two expressions 

Jpdh andJV^/"* dh+C 

Some work will be saved if the first expression can itself be represented by an empirical equa- 
tion directly without having first to determine p and then integrate that result. 
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3. Physically J]*P^ is the amount of air between the earth and the height h in a column 

of unit cross section. This amount can fortunately be read off immediately from a table of 

barometric pressures p at different levels — in fact J*pcEh is simply p a -p 1} the drop in pressure in 
pounds per square foot. 

If F (h) is the amount of air in inches of mercury, the following table gives the value ofF 
from h = 0 to ft = 24,000 feet. 

h =0 . 4,000 8,000 * 12,000 16,000 20,000 24,000 

f: r Q -F^ >*, F=0 4.07 7.68 '.' X0.85 13.62 16.02. 18.00 

.. *. S Now, as the exponential of F (multiplied by a constant) jnust Jbe integrated, it is advisable 

'» i r ;- : ' to have a logarithmic expression to give the empirical relation between h and F. Try: 

TO = 51og 4 (l+aA). 

This holds for fi, = o and has two disposable constants a and &. A least-squares solution could 
be. made to determine the best values for a and b ; but a sufficiently good result may be had 
by passing the curve through or near two sets of values, in the table: 

F, = b log (1 +ah 1 ), F t = h]og(l+a\): 
F, :F 2 =log (l+o \) : log (1 +o A,),l +o \ = (1 +o J,)* 

Let ^=8,000. Then F 3 =7.Q8. Let ^=2^, = 15.36. Then 7i t is just under 20,000 and 
may be taken as ^ = 19,000 by interpolation. Hence" 

1 + 19,000 d = \\ +8000 o)* = 1 + 16,000 0+64,000,000 o» 

and 

0 = 3/64,000. Whence 6=24.1. 

Thus 

FQi) -24.1 log, (1 +37i/64,000) =55.5 Iog„ (1 +3^/64,000). 
The values determined by this function F are 

A = 0 4,000 8,000 12,000 167000 20,000 24,000 
F=0 4.05 7.68 10.8 13.5 15.9 18.1 

These check with the given values to less than 1 per cent, and therefore seem quite good 
enough. 

4. The empirical formula just obtained must be transformed over to suitable units for 
insertion in the integral of the differential equation. Let U be the terminal or limiting velocity 
of the projectile in air of the density p„ at the earth's surface. Then 

cgp Q 17*= Wg or c— W/ U*p 9 
If the units be the pound and foot 



r 



,pShmm ^w x 62,5 x F( ® = 1706 bg « (1 + 3A / 64 > 000 >- 



Taking 0-32.17 and p 0 = . 07608, the result is 

-2g C h jt -1,430,000 , A , 81 \ 

With an error of less than 1 per cent this may be written as 
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An error of 1 per cent seems quite within the range of accuracy possible in such work as this, 
where the variation of atmospheric conditions from the so-called standard table may at any 
. particular time be considerable. 

The integral is therefore 

/1200V /-1200\t 

where the only parameter remaining is the terminal velocity Z7 in air of standard density. 
If in particular a particle starts from rest at the altitude H, 

/1S00V /1300\« 

The integral may be evaluated readily, as the integrand is a simple negative power, but leads 
to a complicated literal formula. 

A simple case would be where U= 1,200. Then 

. _ /. , 3ft \ 64,000 . 



This is the only case in which a logarithm enters; the others are algebraic, e. g., if U= 600, ' 

. , 3ft V 64,000 ^, , 3ft V 8 A , 3ff \~n 

^ = 2 ^V i + 64,000/ 3X3 LA 64,000^ V + 64,000/ J 

5. The formulas thus obtained should not be applied to calculate velocities which are large 
compared with that of sound. It is generally admitted that the simple square law of resistance 
does not hold for velocities much in excess of 800 f t./sec. Just how well the law holds below 
that figure may still be considered doubtful; but two recent authorities (E. Vallier, Balistique 
Exterieure, Encyclopedic des Sciences llathematiques, tome IV, Vol. 6, fasc. 1, p. 15, Gauthier- 
Villars, Paris, 1913; and JL Prescotf, London, Phil. Mag., ser. 6, vol. xx, p. 332, Oct. 1917) 
seem to feel tolerably certain of the law for compact shell-like bodies up to 800 ft./sec. It was 
seen-at the start that irrespective of air resistance, such a velocity would not be obtained from 
rest in a fall of 10,000 feet. If U= 1,200, the velocity of 800 ft./sec. will be reached at the earth's 
surface only when the fall is through some 15,000 feet. (Whether high velocities might not be 
attained at intermediate levels is a question that should not be overlooked.) 

6. It is often stated, as is indeed obvious, that a body falling in a medium of increasing 
density such as the air may reach a maximum velocity and be subsequently retarded before 
striking the earth. The terminal velocity toward which the body strives is of course greater at 
higher levels. The question of reaching a maximum velocity is therefore a question of balance 
between the height of fall and the natural resistance of the body relative to its weight. A 
body of high terminal velocity (in standard air) must fall from a very great height in order to 
attain a maximum greater than the speed with which it reaches the earth, whereas a body of low 
terminal velocity need not fall so far. 

The maximum velocity of a nody falling from rest is obtained by differentiating 

The derivative of an integral with respect to the lower limit is the negative of the integrand. 



1 + 



3g ' 
64,000 



3ft 



I 1+ 64,000 J 
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Hence for the maximum of v, or of v 1 



/1200N* ( 



1 + 64, 000, 



/laoov 

zh \-\nr) 




laooy+r 



-1 + 




= 0 



For a given H this equation will determine tha level h where the maximum velocity is attained. 
If in particular this maximum should be at ~h = 0, the altitude H of fall must be 



For example, if Z7=600, 5=12,500. The maximum velocity will be reached before striking 
the earth only if H> 12,500. If U= 1,200, the bracket has the value e-L= 1.73, and tho.max- 
imum value will not be attained unless the drop is from over 35,000 feet — levels to which the 
empirical formulas used do not remain valid. 

In any case of resisted fall in which a maximum of velocity is attained the value of that 
maximum velocity must- be the value of the terminal velocity at that level, because the body 
has been gaining velocity in the rarer air above and will be losing velocity in the denser air 
below. 

If, therefore, the velocity at &=0 is below the terminal velocity U, it may safely be 
assumed that the value for 7i = 0 has not been exceeded; but if the value for h=0 is larger 
than 27, there has been a still larger velocity at some point of the path. The caution paren- 
thetically suggested at the end of the last article can therefore be observed easily. 

Massachusetts Institute op Technology, September 6, 1919. 
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